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Energies for excited light baryons are computed in quenched QCD with a pion mass of 490 MeV. 
Operators used in the simulations include local operators and the simplest nonlocal operators that 
have nontrivial orbital structures. All operators are designed with the use of Clebsch-Gordan 
coefficients of the octahedral group so that they transform irreducibly under the group rotations. 
Matrices of correlation functions are computed for each UTeducible representation, and then the 
variational method is applied to separate mass eigenstates. We obtained 17 states for isospin 1/2 
and 11 states for isospin 3/2 in various spin-parity channels including = 5/2*. The pattern 
of the lowest-lying energies from each irrep is discussed. We use anisotropic lattices of volume 
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1. Introduction 

The theoretical determination of the baryon spectrum is an important step for lattice QCD. 
Lattice QCD results have provided reasonable agreement with the physical ground state masses of 
different baryons (for a review, see [1]). For excited states of baryons, several groups have reported 
results [2-7], but further improvements are needed for more precise determinations of the masses 
and spins for excited particles. In this report, we show preliminary results for nucleon and delta 
baryon masses including 7^ = 5/2* using improved lattice baryon fields with quenched anisotropic 
lattices. 

Spins of states on a cubic lattice are deduced from the level orderings of energies computed 
in different irreducible representations (irreps) of the octahedral group. Source and sink operators 
have to be designed so as to transform as a definite irrep of the octahedral group in order to identify a 
total angular momentum of a continuum state. There are six double- valued irreps of the octahedral 
group: G\g,Hg,G2g,G\u,Hu and G2„, in which "g" subscript denotes positive parity (gerade) and 
"u" subscript denotes negative parity (ungerade). Subduction pattern of total angular momentum 
to octahedral irreps is summarized in Table |I| 

In general, a lattice operator with a definite irrep couples to multiple states corresponding to 
different total angular momenta according to the subduction pattern. To separate the individual 
energy eigenstates included in an irrep, we employ the variational technique, which involves diag- 
onalization of a matrix of correlation functions. Substantial numbers of basis operators are needed 
in each channel in order to provide sufficiently large degrees of freedom for the determination of 
the spectrum of excited states. Inclusion of multi-hadron operators in the construction of matrices 
of correlation functions is needed to separate the resonances from the hadron scattering states. In 
this work, we only used three-quark operators. Thus the separation of discrete scattering states 
and single-hadron states in the excited spectrum is beyond scope of this work, and we refer to our 
lattice results as "energies" rather than as "masses". 

Table 1: The number of occurrences of double-valued irrep A = {Gi,H,G2} of the octahedral group for 
different values of continuum / (up to 11/2). 
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2. Improved baryon operators 

We construct operators that transform according to each irrep of the octahedral group by em- 
ploying the analytical method based on appropriate Clebsch-Gordan coefficients for the group [8]. 
An alternative method of constructing vast numbers of operators in an automated procedure is pre- 
sented in Ref. [9]. We used complete sets of quasi-local and one-link three-quark operators in this 
work. 
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Quasi-local baryon operators consist of smeared quark fields ^(x, t) that are located at the same 
space point at which the color singlet contraction is taken, 



(IS) 



(IS) (A,X,m)-^,a 
ahc ^a^y 



-dc,c , 



(2.1) 



where i,j,k are color indices, a,b,c are flavor indices, and o;,j3,7 are Dirac indices. Octahedral 
group irreps are denoted by A = {Gig,Hg,G2g,Giu,Hu,G2u}- The corresponding dimensions 
are 2, 4, 2, 2, 4, 2. The rows of an irrep are distinguished by row labels A = 1, • • • ,dA- When 
there are more than one operators in a single (A, A) channel, we use the additional index m to 
label the different operators. These labels are the indices of matrices of correlation functions when 
source and sink operators from a given irrep are used. Isospin I and strangeness 5 are selected by 
coefficients /j^f^- We use the mass for up and down quarks, and do not consider strange quarks in 
this work. The explicit Clebsch-Gordan coefficients c^^''"'' are given in Ref. [8]. 

A simple extension of quasi-local operators is to displace one of the quark fields along a spatial 
axis and to include a straight gauge-link so as to keep the gauge covariance. The displacement 
operator di acted on three-quark operator is defined as follows. 



(2.2) 



where / takes six spatial directions, ±x, ±y, ±z. Assuming that lattice gauge links recover the cubic 
symmetry with large statistics, displacement operators transform amongst themselves under lattice 
rotations. By having the following similarity transformed bases. 
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(2.3) 



the displacements A\,D±fl, and £'0.2 transform according to A|, Ti, and E irreps of the octahedral 
group, respectively. Here Ai is the cubically symmetric, one-dimensional irrep, Ti is the three- 
dimensional irrep that corresponds to £=l and E is the two-dimensional irrep that corresponds to 
^ = 2 for the lowest orbital angular momentum £. These new bases are chosen so that their trans- 
formations resemble those of the spherical harmonics i.e., Ai ~ Yqq, ^+0,- ~ Yn,Yio,Yi-i, 
and £0,2 ^ Y20, {Y22 + Y2-2)- The Ai, Ti, and E one-link operators are now defined as 
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where C I 



A A' A" \ 

X X' X" I ^^^^ Clebsch-Gordan coefficients that takes the direct product of irreps A' 
and A" to yield overall irrep A [8]. 
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We then use improved operators to construct matrices of correlation functions 

C^Kt)8;^'hx'=C^a^;tc%';fl^ma^ (2.5) 

X 

where a displacement operator J, applied to a baryon field denotes a one-link baryon operator for 
/ = zb^;, ±.y, ibz and a quasi-local baryon for / = 0. Lattice operators belonging to different irreps 
or to different rows of one irrep are orthogonal because of the octahedral symmetry of the lattice. 
Different embeddings of a given irrep and row provide operators appropriate to forming a matrix of 
correlation functions, to which the variational method is applied in order to determine the spectrum 



of states. Three Dirac 74 matrices are included in Eq. (2.5) in order to produce a Hermitian matrix 
of correlation functions. 

3. Computational methods 

3.1 Variational method 

Analyses of excited state energies are based upon the matrices of correlation functions of 
Eq. In order to extract the spectrum of energies from the matrix of correlation functions, we 
first average over rows because they provide equivalent results owing to octahedral symmetry. We 
then solve the following generalized eigenvalue equation, 

Lc^£{t)vf{tM = aW(^?o)Icl^)(?o)v?(?,?o), (3.1) 

k' 

where superscript n labels the eigenstates. The symbol C*^^' indicates that the appropriate average 



over rows has been performed. The reference time ?o in Eq. (|3JJ) is taken near the source time 
? = in order to have significant contributions from excited states. The generalized eigenvalues 
a'")(f,fo) are related to the energy by [10] 

a'^"\t,tQ)^e-^"^'-"'\\ + 0{e-\^^\% (3.2) 

where 8E is the difference between and the next closest energy. We have determined effective 



energies En by fitting the generalized eigenvalues to the leading term of Eq. (3.2). 
3.2 Lattice action 

We employed anisotropic lattices with volume 24^ x 64, where the temporal lattice spacing 
at is three times finer than the spatial lattice spacing as. Gauge configurations are generated using 
the quenched, anisotropic, unimproved Wilson gauge action [11] with j3 = 6.1 and the pion mass 
is 490 MeV. A total of 167 configurations is used, but statistical ensembles are effectively doubled 
by use of CPT symmetry. The scale is a^"^ = 6.05 Gev in these lattice, provided by the analysis of 
string tension. 

4. Results and discussions 
4.1 Effective energies 

We have extracted energies for isospin 1/2 and 3/2 channels by diagonalizing matrices of 
correlation functions formed from three-quark operators that share the same octahedral symmetry. 
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We diagonalize matrices of different dimensions by adding or subtracting trial operators to find 
optimal sets of operators. Figure [I] shows low-lying effective energies as a function of matrix 
dimensions. Energies are determined from single-exponential fits to the generalized eigenvalues. 
We found that low-lying energy states are fairly stable with respect to the number of operators 
used, as long as a few important operators are included. It is expected that use of a larger number of 
operators should improves the low-lying spectrum because contaminations from higher-lying states 
are reduced. This expectation appears to hold except for a few excited states that are sensitive to 
the types of operators used. The lowest-lying effective energies in the Gig,Hg,G2g channels with 
isospin 1/2 are plotted in Fig. 0. 
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Figure 1: Effective masses in lattice units vs. dimension of matrices of correlation functions. Upper panels 
show the positive-parity channels and lower panels show the negative-parities channels, no matrices of 
correlation functions are obtained for the A,G2g/„ irreps. 
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Figure 2: Effective energies for the lowest states in Gig^H^ and irreps with isospin 1/2. Time ranges 
used for and energies obtained from fits are shown by horizontal lines. 
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4.2 Pattern of low-lying energies 

Figure |3| show the lowest energies from each irrep in isospin 1/2 and 3/2. Energies are con- 
verted to physical units. They are determined based the diagonalization of matrices of correlation 
functions consisting of the best sets of operators. This graph shows the pattern of level orderings for 
the different irreps, but it is not to compare the results with experimental values because our pion 
mass is too heavy and there may be significant discretization errors. Nevertheless, it is interesting 
to compare the pattern of level orderings with empirical data. 

In the positive-parity channels, besides the nucleon ground state and the delta ground state, 
isospin 1/2, Gig and Hg states are nearly degenerate within errors. If this degeneracy holds in the 
continuum limit and no degenerate partner exists in G\g, these two eigenstates correspond to the 
subduction of = 5/2+. In our simulations we have no evidence for a 3/2+ state that has lower 
energy than 5/2+ state. Indeed in nature N{5/2'^, 1675) is lighter than A^(3/2+, 1705). Our results 
agree with this pattern, though our simulation conditions are far from physical. Previous lattice 
calculations that used Hg{Hu) operators, such as Rarita-Schwinger projected operators, assumed 
that the lightest state corresponds to spin 3/2+ (3/2^). Spin identification of the lowest H state is 
uncertain unless G2 irreps are included, i.e., the state could have spin 5/2 or higher. The G2 baryon 
operators cannot be constructed by local operators, therefore displaced operators are needed for the 
lattice simulations of excited baryon spectrum. 
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Figure 3: The energies obtained for each symmetry channel of isospin 1/2 and 3/2 baryons are shown. The 
energies are converted in MeV using the scale a,^' — 6.05 GeV. 
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In the negative-parity channels, the level orderings from the lattice simulations go as follows: 
the N*,Hu corresponding to 7^ = 3/2^ (no degenerate partner with N*,G2u) is the lowest, the 
N* ,G\u corresponding to 1/2^ is slightly above, the A*,Gi„ corresponding to 1/2^ is next, the 
A*, Hi, corresponding to 3/2^ (no degenerate partner with A*,G2u) is slightly above, the N*,G2u 
corresponding to 5/2~ is next and finally the A*,G2u corresponding to 5/2^ is the highest. These 
spin identifications are only tentative. In our simulations, we cannot find more that one operator in 
isospin 3/2, G2g/u chaimel, thus there is uncertainty in these chaimels. The pattern of level orderings 
in the negative-parity charmels is consistent with the experimental data, except for the ordering of 
the A*,Hu and the N*,G2i, states. 

It is intriguing to note that comparison of energies between different parities gives very differ- 
ent conclusions. For instance, the (lowest) A*, 3/2^ state should appear between the N*,5/2'^ and 
the A*, 1 /2"'", but our calculation does not agree with this pattern. We do not find a positive-parity 
excited state that has Ughter energy than the lowest negative-parity state in isospin 1/2. Further 
studies of lattice simulations are needed to understand realistic mass gaps between parities. 

For more details about the quenched baryon spectrum based on extended operators, readers 
are referred to Ref. [12]. 

This work was supported by the U.S. National Science Foundation under Award PHY-0354982 
and by the U.S. Department of Energy under contracts DE-AC05-06OR23177 and DE-FG02- 
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